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SUMMARY 


Some measurements were made of the aerodynamic forces and moments 
acting on a rectangular wing of aspect ratio 2 which was oscillated about 
the midchord. These measurements were made at four frequencies (31, 43, 

54, and 62 cps) over a range of Mach number from O.I 5 to 0.8l, a range of 
reduced frequency from 0.15 to 1.32, and a range of Reynolds number from 

0.60 X 10^ to 9.21 x 10 6 . It was feasible to compare results of a por- 
tion of these measurements with some published experimental data and, in 
general, reasonable agreement was found to exist. 

An appendix was presented to show the correction of the root reaction 
for inertia and aerodynamic effects in order to determine the total aero- 
dynamic load. 

. Comparison of the measured aerodynamic forces and moments with those 
predicted by the method of Reissner and by the method of Lawrence and 
Gerber for wings of aspect ratio 2 in incompressible flow showed generally 
good agreement. Comparison of the measured quantities with those pre- 
dicted by two-dimensional-flow theory indicated that the effects of finite- 
ness of span on the aerodynamic forces and moments are considerable . 

Some experimental results pertaining to the influence of wind-tunnel- 
wall effects on nonsteady aerodynamic measurement have been included in 
an appendix. 


INTRODUCTION 


A need exists for experimental measurements of oscillating air forces 
because of the significance of these forces in flutter and related fields, 
and in order to assess past and present theoretical work. Despite the 
importance of the problem, only a limited amount of data, both experimental 
and theoretical, exists for restricted ranges of aspect ratio, Mach number, 
and Reynolds number. (See, for .sl-nSibanc e . ref. 1.) 
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Theoretical treatment of aspect-ratio effects on the oscillatory 
aerodynamic coefficients for subsonic speeds is in a state of develop- 
ment and is not yet in a form convenient for numerical comparison. There 
is difficulty in mathematically representing the physical phenomena and 
certain assumptions necessary to obtain a solution are doubtful, particu- 
larly those associated with tip effects. Consequently, current experi- 
mental measurements will be compared with readily available experimental 
and theoretical oscillatory coefficients . 

This paper presents some experimental measurements of the oscilla- 
tory aerodynamic forces and moments acting on a rectangular wing of 
aspect ratio 2 which was oscillated about the midchord. These coef- 
ficients are presented for a range of reduced frequency from 0.15 to 
1.32 and for a range of Mach number from 0.15 to 0.812. The Reynolds 

number ranged from 0.60 x 10^ to 9-21 x 10^. These measurements were 
made by using a resonant oscillation technique in the Langley 2- by 
4-foot flutter research tunnel with air or Freon-12 as a testing medium. 
A comparison of the measured values has been made with some existing 
experimental data and with the results of the analyses of references 2 
and 3, that is, of Reissner and of Lawrence and Gerber. In order to 
establish some convenient reference values, coefficients for two- 
dimensional incompressible flow (ref. 1) are indicated. The results 
are presented in tabular form for quantitative evaluation and in graphi- 
cal form for qualitative examination and comparison. 


SYMBOLS 


A 

a 

s 

c 

St 

Ss 

H 

k 


aspect ratio 
speed of sound, fps 
semispan of wing, ft 
chord of wing, ft 

logarithmic damping coefficient of wing in airstream 

logarithmic damping coefficient of wing in a near vacuum 
tunnel height, ft 

reduced-frequency parameter, u>c/2v 
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K g effective spring constant of oscillating system, ft-lb/radian 

L oscillating lift vector acting on the wing for oscillations 

about the midchord axis, positive when acting downward, 

L = |L|e i ( t0 " trK P) = -(npv^sc|a|)(li + il2) eiayt 

nondimens ional coefficient of lift in phase with the angular 
displacement 

Z ‘2 nondimens ional coefficient of lift in phase with the angular 

velocity 

<f> phase angle between lift vector and incidence vector, 

tan -1 (Z 2 /Zi) , deg 

Mo, oscillating moment vector acting on wing for oscillation about 

midchord, referred to the axis of rotation, positive for 
leading edge up, 

Mo, = Jl^j e^ 03 "*^ 6 ) = («/2)pv^sc^ |a| (mi + 

mi nondimensional moment coefficient in phase with the angular 

displacement 

I s effective inertia of the oscillating system accounting for 

the dynamic deformation of the system, ft-lb-sec2/radian 

m 2 nondimensional moment coefficient in phase with the angular 

velocity 

9 phase angle between moment vector and incidence vector, 

tan"-*-(m 2 /mi) , deg 

M Mach number 

R Reynolds number 

v velocity of test medium, fps 

a angle of incidence at s/2 span station as a function of 

time , | a | e , radians 

p mass density of test medium, slugs/cu ft 

circular frequency of oscillation of wing, radians/sec 



o>i 
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co^ circular frequency of first natural wing bending, radians/sec 

o> vac circular frequency of oscillation of wing in a near vacuum, 

radians / sec 


APPARATUS AND METHOD 
Tunnel 


The Langley 2- by 4-foot flutter research tunnel was used for the 
tests reported herein with the test section modified to be rectangular 
in shape, measuring 4-5.75 by 24 inches. The test mediums were air and 
mixtures of air and Freon-12 as noted for each set of data. The use 
of the air and Freon-12 mixture as a test medium permits the attainment 
of approximately twice the reduced frequency obtained in air for a given 
Mach number and frequency. The choking Mach number for these tests was 
approximately 0.92. The wing was mounted in the test section as shown 
in figure 1 and with its wall reflection had an aspect ratio of 2 . 


Wing Model 

The semispan wing -model had a rectangular plan form with a 12-inch 
chord and a 12 -inch semispan corresponding to an aspect ratio of 2 . 
Fabricated construction was employed; a steel box spar carried four 
evenly spaced ribs to which plywood skin was attached forming an NACA 
65 AOIO airfoil section. The wing was designed to have high natural 
frequencies in order to reduce the amount of correction to the measured 
forces due to elastic deformations and to bending inertia loads . The 
first natural cantilever bending frequency ranged from 125 to 130 cps . 

The semispan wing model was mounted as a cantilever beam at the 
tunnel wall in an oscillator mechanism. The mount permitted the wing 
to oscillate in pitch about the midchord axis . The wing was mass balanced 
about this axis of oscillation in such a way that there were no lift 
reactions when the wing was oscillated in a near vacuum. 


Oscillating Mechanism 

The oscillating mechanism may be considered as a simple torsional 
vibratory system as illustrated in figure 1. The system consists of a 
torsion spring which is fixed at one end, a hollow steel shaft which is 
supported by bearings, and the semispan wing which has a base plate that 
is flush with the tunnel wall. The mechanism was oscillated in torsion 
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at the natural frequency of the simple spring-inertia system by applying 
a harmonically varying torque through the shaker coils attached to the 
shaft. Different springs were used to permit a choice of frequency of 
oscillation. The springs used and the resulting torsional natural fre- 
quencies in a vacuum were as follows: 


Spring 

°Vac> 

radians /sec 

A 

31 X 2rt 

B 

43 

C 

54 

D 

62 


The bearings were contained in housings which were carried on 
columns. These columns were designed to include stress -sensitive regions 
and were equipped with strain gages from which the aerodynamic lift could 
be determined. The vertical reactions at the fixed end of the torsion 
spring were negligible because of the rigidity of the steel tube and the 
small deformations experienced by the strain-gage columns. 

The electromagnetic shakers consisted of stationary coils furnishing 
a steady magnetic field and moving coils which were attached to the steel 
shaft. The moving coils were driven by a variable-frequency oscillator. 
The moving coils were alined so that the direction of the applied force 
was perpendicular to the direction of the lift. Provision was made for 
interrupting the power to the moving coils in order to obtain a power-off 
decaying oscillation. 


Instrumentation 

The instrumentation of the experiment was designed to provide signals 
that were a measure of the lift and angular position at any instant and 
to provide a means of measuring their amplitude and time relationship. 

The lift reactions were converted to electrical signals by means of wire 
strain gages attached to the supporting columns. The gages were connected 
so that only lifting loads were sensed. An electrical signal from a wire 
strain gage mounted on the torsion spring so as to sense torsional strains 
was calibrated to give the angular displacement in terms of the wing inci- 
dence. The signals were filtered and measured with vacuum-tube voltmeters 
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The angular-position signal was recorded on a recording oscillograph 
during the decay of the oscillation for the purpose of obtaining the 
damping factor. 

The phase measurements were made with an electronic counter chrono- 
graph. The time lapse between a given point on a lift signal and a 
corresponding point on the position signal was measured while the wing 
was oscillating at a constant frequency. The period of oscillation was 
measured with this instrument by determining the time lapse between 
corresponding points on the same signal. 


Calibration 

/ 

The angular position of the wing was dynamically calibrated with the 
signal from the torsion strain gages by a photographic technique. Time 
exposures were taken of a fine chordwise line on the tip of the wing for 
various amplitudes while the signal output was recorded. The amplitude 
of oscillation of the wing was obtained from the enveloped position of 
the line on the wing tip and correlated with the strain-gage signal. By 
using this procedure and a line on the leading edge, it was determined 
that, at the maximum frequency of oscillation (62 cps), the tip angle of 
incidence exceeded the root single of incidence by less than 2 percent. 

The signals from the balsmce columns were calibrated in terms of 
pounds of force per unit of signal strength. Known loads were applied 
to the wing, and the column reactions were determined by treating the 
wing shaft system as a simple beam with overhang. (See, for instance, 
fig. 1.) 

The reaction forces were then related to the respective signals . 

The meters were calibrated dynamically by using a voltage divider referred 
to the open-circuit calibration of the strain gages. The vacuum-tube - 
voltmeter readings are believed to be within ±4 percent of true signal. 

In order to minimize errors in phase <p introduced in the electrical 
operations, a tare value of phase was obtained at each reading by applying 
either the lift or incidence signals through both channels of the elec- 
trical circuits. The phase -measuring system was calibrated at various 
frequencies by using standard resistance-capacitance phase-shift circuits, 
and by using a cam-operated set of cantilever beams on which strain gages 
had been mounted. The latter system had distortion and noise, approxi- 
mating the worst tunnel condition. Calibrations of the phase meter indi- 
cated that the phase angle may be determined within -3° of true value 
with a noisy signal and within *0-5° of true value with a clean signal. 
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Data Reduction 

The lift forces as received from the balances contain an aerodynamic 
component and an inertia component which arises from the bending deforma- 
tion of the wing. In order to correct the measured lift to the aero- 
dynamic lift, it was necessary to correct for the inertia forces due to 
wing deformation. A discussion of this correction is given in appendix A. 
The inclusion of this correction leads to a factor which, when multiplied 
by the measured lift, gives the actual applied lift. Values of the factor 
are O .98 for spring A, 0.95 for spring B, 0.91 for spring C, and O .87 
for spring D. In order to estimate possible error incurred by neglecting 
the aerodynamic forces and moments arising from the bending deformation, 
these forces were included in the analysis in appendix A and were found 
to be less than 1 percent of the correction due to wing deformation caused 
by the inertia forces. The phase angle 9 contained a component due to 
these forces that tended to increase cp by less than 1°. The moments due 
to the bending deflection were found to be negligible relative to the magni- 
tude of the measured moment. As the aerodynamic effects due to wing bending 
were within the accuracy of the measurement, no effort was made to adjust 
for these quantities. 

The inphase moment was determined from the change in resonant fre- 
quency due to air flowing over the wing as indicated in reference 5 * 

Since the torsional damping was relatively small, its effect on the fre- 
quency is neglected and the entire shift is attributed to the inphase 
moment. The coefficient of the inphase moment is given by 


mi 


-K c 



rt 

2 


pv 2 


c^s 


\^vac J 



The dependency of m^ on the small difference of two quantities of the 

same magnitude leads to considerable loss in accuracy and consequent 
scatter in the data. This scatter became so large for the torsional 
springs C and D that those coefficients are not presented. 

The quadrature -moment coefficient was determined by operating on 
the logarithmic decrement of the power-off decay of the oscillation. The 
coefficient was given by 


-Hl2 = 


w 

— pv 2 C 2 S 
2 


et - 



The derivation of this equation is treated in appendix B. 
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The phase angle 0 between the moment vector and the angle of inci- 
dence was obtained from the relationship 9 = tan -1 (m2^m-]J . The lack of 
precision of determining m^ and mr> directly affects the degree of 

accuracy of 0; the values of 0 are not expected to be more accurate 
than its components. 


RESULTS AND DISCUSSION 


The experimental results for the measured aerodynamic forces and 
moments are presented in tables I to IV. These results cover four over- 
lapping ranges of reduced frequency (one range in each table) because 
the Mach number was varied over a range while the frequency varied only 
a small amount (due to the change in aerodynamic moment ) . The measured 
values of lift and moment coefficients and their respective phase angles 
are shown in tables I and II, whereas the inphase moments and moment phase 
angles are omitted from tables III and IV because of the lack of these data. 

The calculated values of the various theories are given in tables V 
to VII for convenience of discussion and comparison with experimental 
coefficients. Some calculated values are also shown in various figures 
for ease in comparing trends. Experimental data relating to the influence 
of wind-tunnel walls on the aerodynamic forces and moments are presented 
in table VIII. 

Before presenting the actual coefficients of lift and moment and 
their respective phase angles, the effect of Mach number and Reynolds 
number on the coefficients will be briefly mentioned. A comparison of 
some of the current data will then be made with some other experimentally 
determined coefficients obtained from another source (ref. 6). Compari- 
sons of the current data with the two aspect-ratio theories will then be 
made. For reference, coefficients for two-dimensional incompressible flow 
will also be shown. The section will be concluded with a brief discussion 
of wind-tunnel-wall effects. 

A study was made to determine the effects of Mach number and Reynolds 
number on the aerodynamic coefficients. Since the testing technique used 
did not readily permit holding either k, M, or R constant and varying 
the remaining two parameters, a considerable amount of cross plotting was 
necessary to obtain any indication of an effect due to Mach number or 
Reynolds number. It was found that for this purpose there were insuffi- 
cient data and much extrapolation and interpolation were necessary. In 
light of the experimental inaccuracies (perhaps of the order of the par- 
ticular effects sought) and the operations necessary to obtain the results, 
no quantitative information could be obtained. This study indicated that 
the aerodynamic coefficient possibly was influenced to some extent by 
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both Reynolds number and Mach number. However, for the ranges of speed 
and frequency covered in this series of experiments, the over-all effects 
of the Reynolds number and Mach number appear not to be of first order, 
and to lie perhaps within the accuracy of the experiment. 

In view of the paucity of experimentally determined oscillatory aero- 
dynamic coefficients for finite wings, it is of interest to compare these 
data with available data from other sources. A recent paper, reference 6, 
has presented experimental data for a rectangular wing having the same 
aspect ratio and axis of rotation as the wing presented herein, but for 
a lower range of Mach number and Reynolds number. It should be r em arked 
that there are differences in techniques; the moment coefficients of 
reference 6 were obtained for the case of steady oscillation, whereas 
the moment coefficients of the present paper were obtained from the 
method of decaying oscillations. However, in reference 5 it was shown 
that for a two-dimensional wing the damping moment obtained from steady 
oscillations was in agreement with that obtained from decaying oscillations. 

A comparison of the results presented herein with those published in 
reference 6 is shown in figure 2. It may be seen that, for the magnitudes 
of the lift and moment, there is good agreement. This agreement is grati- 
fying and, since different methods of measurement were used, the magni- 
tudes shown may be considered valid. Large discrepancies may be noted 
between moment phase angles. The data of reference 6 fall close to a 
phase angle of 0^ and show some values which indicate negative aerodynamic 
damping. The reasons for the discrepancy in the moment-phase-angle 
results between reference 6 and the present paper are not known. A com- 
parison of current moment-phase angles with results of calculations will 
be given later in this paper. 

Additional data on the experimental damping-moment coefficients have 
been given in reference 7. These coefficients were obtained by the method 
of decaying oscillations for a rectangular wing of aspect ratio 2 with 
pitch axis at the midchord. The results of reference 7 are compared with 
the data of this paper in figure 3, the solid and dashed curves being 
taken from reference 7 and the symbols representing results from this 
paper. Despite the differences in Reynolds number, it may be noted that 
the data from the present investigation are in basic agreement with other 
published results. 

The comparisons of the experimental coefficients and phase angles 
with the theoretical oscillatory coefficients and phase angles for a 
finite wing will be made only for the incompressible case in this paper. 

The results of two theories, namely, those of Reissner (ref. 2) and 
Lawrence and Gerber (ref. 3), were readily available, whereas coeffi- 
cients for finite wings in compressible flow were not. In making these 
comparisons, it is recognized that Reissner limits the applicability of 
his theory to aspect ratio 3* No effort has been made in this paper to 
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evaluate the relative merits of the two theories. In order to illustrate 
the influence of finiteness of span, further comparisons of the experi- 
mental data have been made with the coefficients for a two-dimensional 
incompressible fluid (ref. 4 ) . These comparisons may he made by referring 
to figures 4 to 8. 

In figure k may be seen the coefficient of the magnitude of the lift 
vector as a function of l/k. In this presentation it should be recalled 
that the data for different torsional springs at a given value of l/k 
were obtained at different Mach numbers. It may be seen that, for the 
low values of l/k, experimental data fall between the theory of Lawrence 
and Gerber and the theory of Reissner, whereas for the higher values of 
l/k, the theory of Reissner appears to be in agreement with the data and 
the theory of Lawrence and Gerber falls below the measured lift coeffi- 
cients. The higher M values correspond to the higher l/k values. 

The steady-state lift coefficient corrected for aspect ratio by the 

A A 

factor for the incompressible case and by the factor 

A + 2 J i — 

A Ml - M2 + 2 

as suggested in reference 8 for the compressible case evaluated for 
M = 0.70 also are shown in figure 4. The value of the incompressible 
approximation to the steady-state lift coefficient is a fair representa- 
tion of the experimental values for l/k greater than 1.5. The experi- 
mental values do not appear to vary appreciably from the value determined 

A 

by the correction . Thus, it might appear that the steady-state 

A + 2 

lift coefficient could be used as a basis for est ima ting the magnitude 
of the lift coefficient over a considerable range of reduced frequency. 

The inadequacy of estimating these coefficients with two-dimensional 
incompressible-flow theory with no correction or modification is also 
indicated in figure 5* This correction would appear to vary for this 
case from O.67 for a value of l/k of 2.5 to O.56 for a value of l/k 
of 6.0. 

The measured phase angle of the lift Vector for spring A is com- 
pared with theoretical phase angles in figure 5. The theories may be 
noted to agree in general and the theory of Reissner is in good agree- 
ment with the measured values, the theory of Lawrence and Gerber falls 
above the measured values by a few degrees, and the two-dimensional- 
flow theory is seen to fall below the measured values by a few degrees. 

The phase angle as given by the aspect-ratio theories and as measured 
appears to approximate a linear function of reduced frequency in the 
range of k from 0.2 to 0.7» Since the frequency was essentially 
constant for a given torsional spring, the low values of reduced fre- 
quency were obtained at the higher values of Mach number. In general, 
the phase angle of the lift appears to be predicted by the aspect-ratio 
theories of references 2 and 3- 
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Although some of the moment coefficients have already been pre- 
sented as magnitudes and phase angles (fig. 2), the components of the 
measured moments are compared with the components of the theoretical 
moments. It is felt that, since the measured values were determined as 
components, it is appropriate to conpare these values with theoretical 
components. By referring to figure 6 it may be seen for the Inphase- 
moment coefficient that the aspect-ratio theories give values which are 
in general agreement with the data. The two-dimensional-flow theory is 
given for reference only. With regard to the values given by the aspect- 
ratio theories, it is seen that there is little difference between them 
although the theory of Lawrence and Gerber appears to agree a little 
better with the measured values. One other comparison that may be made 
is with the steady-state value of the moment coefficient, for which there 

is no damping moment. The coefficients with the correction — - — and 

A + 2 

including compressibility — A are indicated in figure 6. It 

A\J 1 - M 2 + 2 

appears that the steady-state coefficient uncorrected for Mach number 
but corrected for aspect ratio could be used as a basis for estimating 
the magnitude of inphase -moment coefficients over the range of l/k 
greater than 1.5* 

With regard to the damping-moment-component coefficients, these 
data and the theoretical coefficients are shown in figure 7- The two- 
dimensional -flow theory is given again for reference only and, of course, 
needs large correction factors to make it apply to the data. The scatter 
in the data in conjunction with the possible influence of Mach number and 
of Reynolds number precludes exact conclusions with respect to the agree- 
ment of the data with the theoretical coefficients of the aspect-ratio 
theories. Although each of the aspect-ratio theories is to some extent 
in agreement with measured coefficients for the damping moment over some 
range of l/k, it may be noted that neither aspect-ratio theory covers 
the over-all range adequately. 

The phase angle between the moment vector and the angular -position 
vector was obtained from the ratio of the out-of -phase- and inphase- 
moment coefficients and is shown in figure 8. The scatter in the data 
is attributed to the scatter present in the components of the moment 
coefficients. It may be seen that the measured coefficients are in 
fair agreement with the coefficients given by the aspect-ratio theory 
of Reissner for l/k greater than 2, whereas the theory of Lawrence 
and Gerber gives slightly lower values. In this instance, it appears 
that the two -dimensional -flow theory predicts substantially the same 
values as the aspect-ratio theories although it is slightly different 
in trend. 
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Before closing the discussion on the measured aerodynamic coeffi- 
cients, it is appropriate to mention the possible influence of wind-tunnel- 
wall effects on the measured quantities. A recently published analytical 
investigation (ref. 9 ) of the effects of wind-tunnel walls on air forces 
on a two-dimensional oscillating wing at subsonic speed demonstrated the 
possibility, under certain conditions, of the existence of large t unn el - 
wall effects, associated with an acoustic resonance phenomenon. It was 
also pointed out that similar conditions exist for three-dimensional 
flow. It was shown in this reference that a condition for a maximum of 


distortion is satisfied by the equation M cr = 




1 - 


cu 2 H 2 

rt 2 a 2 


vhere 


M CI 


IS 


the Mach number corresponding to the circular frequency of oscillation 
and tunnel height H at which the phenomenon will occur. The symbol a 
is the speed of sound.. At the present time, no quantitative calculations 
for a finite wing are available . 


For the purpose of showing the proximity of the data reported herein 
to the region of critical tunnel-wall interference based on two-dimensional 
flow, plots of k against M for the various torsion springs are shown 
with a curve of critical tunnel-wall effects in figure 9* The curves 
marked A, B, C, and D are well away from the curve of critical wall 
effects and, thus, the tunnel-wall effects might be expected to be small. 
The range of critical M cr for the data given in the body of this paper 

is between O .89 and O. 96 , whereas the highest Mach number reported is 
0.8l. The curve marked D' intersects the critical curve at a value 
of M = O.U 7 or k = 0.77 and represents data from spring D in Freon-12. 
Some evidence of the tunnel-wall effects based on the D' curve is shown 
in appendix C. 


CONCLUDING REMARKS 


Some measurements of the oscillating aerodynamic forces and moments 
acting on a rectangular wing of aspect ratio 2 oscillating about the mid- 
chord were made at four different frequencies. The measurements were 
made over a range of Mach number from 0.15 to 0.8l, a range of Reynolds 

number from 0.60 x 10^ to 9*21 x 10^, and a range of reduced frequency 
from 0.15 to 1 . 32 . 

Appendixes are presented to show the correction of the root reaction 
for inertia and aeroelastic effects in order to determine the total aero- 
dynamic load, to show the determination of the aerodynamic damping -moment 
coefficient, and to show some experimental evidence of wind-tunnel-wall 
effects on an oscillating wing. 
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A study of the effects of Mach number and Reynolds number indicated 
that, in the range of the experiment, the over-all effects appeared to 
be small, although there were insufficient data even to determine quali- 
tatively their trends. The phase angle between the lift vector and angu- 
lar position was seen to vary fairly linearly over the range of reduced- 
frequency parameter studied. The inphase lift and the inphase moment 
remained essentially constant with frequency parameter, whereas the quad- 
rature lift and quadrature moment were found to increase with an increase 
in frequency parameter. 

Comparisons of the data were made with existing published data and 
with theoretical incompressible coefficients obtained from the aspect- 
ratio theory of Reissner, the aspect-ratio theory of Lawrence and Gerber, 
and the two-dimensional-flow theory. A comparison of the experimental 
data of this paper with other experimental data showed that good agree- 
ment was obtained for those coefficients that could be accurately deter- 
mined. In the comparison of these data with the theoretical coefficients 
it was found, as might be expected, that the two -dimensional -flow theory 
was inadequate for predicting the experimental coefficients. However, 
the coefficients given by the aspect-ratio theories were generally in 
good agreement with the experimental data. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 17, 1953* 
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APPENDIX A 

CORRECTION OF ROOT REACTION FOR INERTIA AND AEROELASTIC 
EFFECTS TO DETERMINE THE TOTAL AERODYNAMIC LOAD 


In the present experiment on a cantilever wing, knowledge is desired 
of the total aerodynamic load which develops solely from the torsional 
oscillations of the wing. This aerodynamic load is not equal precisely 
to the reaction at the wing root because of the presence of secondary 
bending reactions which come as a result of the freedom of the wing to 
be excited slightly in a bending oscillation. A correction must, there- 
fore, be applied to the measured root reaction to obtain the aerodynamic 
load associated directly with the torsional motion. This appendix derives 
and shows the magnitude of this correction. The derivation is made in 
general terms of a wing of variable cross section; the correction is then 
applied to the uniform wing. 

On the basis of the engineering beam theory, the differential equa- 
tion for bending of the wing is 


( i + ib ) & ei 3 =p 


(Al) 


where g is the structural damping coefficient and p is the intensity 
of the applied loading. With the choice of a strip analysis approach, 
the loading for the case under consideration may be written 


P = -my — y - | pcv(F + iG)y + P 


(A 2) 


The first term in this expression is the inertia force associated with 
the wing mass; the second and third terms refer, respectively, to the 
"apparent air-mass" inertia effect and the' aerodynamic damping associated 
with bending oscillations; and the fourth term refers to the torsionally 
induced aerodynamic loading, which here is regarded as the "applied forcing 
function." The second and third terms were established by using oscillating 
flow theory for two-dimensional incompressible flow as a guide; the lift- 
curve slope a and the F and G coefficients, which are like the inphase 
and out-of -phase Theodorsen flutter coefficients, are to be selected as 
appropriate to the case being treated. 
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Since harmonic motion is involved, the loading and the deflection 
may he written 

P = L(x)e icut (A3a) 

y = Y(x)e icut (A3b) 

Now, with the use of equations (A2) and (A3), equation (Al) reduces to 

(l + ig)^— El = <r£mY - i §• pcva>(F + iG)Y + L (A4) 

dx 2 dx 2 2 


where 


m 


m + 


Jtpc2 

~ 


A convenient and fairly accurate approximate solution to this equation 
can he obtained by expressing the deflection in terms of the fun dam ental 
bending vibration mode of the wing, the choice here of only a single mode 
expansion being considered adequate since the forcing frequencies used 
in the experiments were below the fundamental wing frequency. Thus, 

Y = a lYl (A5) 

where represents the response amplitude to be determined and y± 

is given in terms of unit tip amplitude and satisfies the equation 


d2 

dx 2 


__ o- 

EI = oq^my^ 

dx 2 


(A6) 


In accordance with the Galerkin procedure for solving differential equa 
tions, equation (A5) is substituted into equation (Ah) which is then 
multiplied by y 1 and integrated over the length of the wing. The 
result, with the use of equation (A 6), is 
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ap (l + ig)a^M-^ = apa/^Mp 


ia-L | ?L u£(F + iG)Ap + J Lyp dx 


(AT) 


where 


% « r my p^dx 
J 0 


Ai = f — yi 2d ^ 

Jq c r 


m^* — 


k = 


Jtpc r £ 

“ T 

<nc r 

2v~ 


and c r is some convenient reference chord, usually taken at about the 
3/^-span station. The desired response amplitude can now be determined 
directly from equation (A7 ) ; hence, 


a l “ * 

(aq_2 - ur)Mp 



£ | + i^SMi + | | urAm?) 


(A8) 


The loading on the beam is now written in terms of aj_. 
tion of equations (A3) and (A5) into equation (A2) gives 


Substitu- 



o^yi - i | pcvo)(F + iG)y 1 ] 


icst 


P = 


a^ + L > e 
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The root shear or reaction, which may be designated V 0 e ia> ^, may be 
found by integrating this loading over the length of the beam} thus, 

A 


e iaJt = „ a^N-i - i(F + iG) - B-. 

° 1 it k x 


& 1 + 


L dx 


'0 


? e 


ia>t 


(A9) 


where 


Nq = / %q dx 

'-'0 


B i = / — yi 

Jq c r 


Substitution of equation (A8) into equation (A9) and cancellation of the 
harmonic terms gives 


v o - 


1 , (°1 + ic 2 ) d 


D 1 + 1D 2 


L dx 


(A10) 


where 


c = Hi + - 2 rc^l 

^ M]_ it k Mq 


a F % B 1 


C2 = - - f 

Jt k Ml 


p. _ (D L a G ^1 

1 ^2 it k Ml 


D2= g l*if^ 


o£ It k Mx 


d - lo ^ 4 Jo L 


(All) 
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As stated earlier, knowledge is desired of the total aerodynamic 
load that is associated with the torsional oscillations of the wing. 
This total load is found directly from equation (AlO) to he 



L dx = 


D 1 + iD 2 
C-j_d + + i (c 2 d + D 2 ) 


This equation, when expressed in the complex notation of a modulus and 
a phase angle, becomes 



L dx = 


»1 2 


(Old + Vi)‘ 


+ Dgg 

+ (C 2 d + D 2 ) 2 




(A12) 


where 


= tan 


-1 


(Cl^2 “ C 2 Di)d 


d (CiD]_ + C 2 D 2 ) + D x 2 + I> 2 2 


(A13) 


This is the final equation sought. Thus, the magnitude of the torsionally 
induced air load is found simply by multiplying K by the magnitude of 
the measured root shear* j the phase angle between this load and the root 
shear is given by A word about the coefficients Cq, 1^, and d 

may now be in order. All terms in these coefficients which contain the 
lift-curve slope a are related to the aerodynamic damp-trig effects 
associated with the bending oscillations. A comparison of the second 
term with the first term in D 2 , for example, will indicate how strong 

the aerodynamic damping is in relation to the structural damping. The 
nondimens ional term d may be seen to depend on the distribution of the 
air load L(x). For most practical cases, it is considered sufficiently 
good to evaluate this factor on the basis that the air load has 
an elliptic distribution. 

Some simplification results when the aforementioned relations are 
applied to the case treated herein, that is, to a uniform cantilever 
wing. In this case, it is easily shown that 
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Bl 


y]_ dx = 0.39157 


J 


Hence , 


ci- 1.566(1 ♦!§!) 


c 2 = -1.566 a 1 2* 

* « k ffl 


9 


D. = — 

1 CD 2 


, a G “r 


(A15a) 


T, , a F “r 

2 8 o£ * k m 


and, for an assumed elliptic loading. 


d = 0.29 


(A15b) 


In order to determine the correction K and the phase angle 0^ 
for the four different spring combinations that were used in the tests. 
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equations (A12) and (A13) were used together with equations (A15). The 
structural damping coefficient used for the wing was g = 0.008. The 
lift-curve slope was taken equal to the theoretical value of 2 n multi- 
plied by the often used aspect-ratio correction — — — , and was thus 

A + 2 

taken as it since the wing is of aspect ratio 2. The F and G func- 
tions were arbitrarily chosen as those for two-dimensional incompressible 
flow. The results obtained for the four cases are shown in the following 
table : 


Torsion 

spring 

co x (2it) -1 

Structural damping 
only, a = 0 

Both structural and 
aerodynamic damping 

K 

01 > deg 

K 

01, deg 

A 

31 

O .98 

0.012 

0.98 

O.06 

B 

43 

•95 

.030 

.95 

.11 

C 

54 

• 91 

• 060 

.91 

• 19 

D 

62 

• 87 

.093 

.87 

•30 


In order to gain an insight as to how aerodynamic damping effects 
compare with structural damping effects, the calculations were made for 
two conditions: (l) with structural damping only and ( 2 ) with both 

structural and aerodynamic damping included. The table shows the results 
for these two conditons. No differences are noted in K for these two 
conditions. However, K decreases from a value of O .98 to a value of 
O .87 as the forcing frequency increases, thereby denoting a correction 
that ranges from 2 to 13 percent. Although a difference in phase angle 
is noted for the two damping conditions, the important thing to note 
is that in all cases the phase angle is a negligible quantity. 
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APPENDIX B 

DETERMINATION OF THE AERODYNAMIC DAMPING-MOMENT COEFFICIENT 


In this appendix the method used in obtaining the aerodynamic 
damping-moment coefficient from the power-off decaying oscillations of a 
torsional spring-inertia system is given. By assuming a linear problem, 
particularly with respect to the aerodynamic coefficients, and by using 
the concept that the structural damping moment is in phase with the 
angular velocity but proportional to the angular displacement, the 
differential equation of motion of the system is given by 

I s a + k s /l + ig s ^a = £ pv^sc^a^ + im. 2 ) (Bl) 


or 


a + 


c + IB 
A 


0 


By definition. 


and 


k sS s - | pv^sc^amj 
A = 5£ 



B 

C 


k sSs - 5 pv 2 c 2 sm 2 

A 

k s “ | pv 2 c 2 sm 1 


Equation (Bl) has a solution of the following form (see, for instance, 
page 86 of ref. 10): 


a = 




2A 


a | e ’ 


- \(l+A 2 ) 1 / 2 -l+i\ ( 1+A 2 ) 1//2 +l 


t 


(B2) 
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For small values of 


A, equation (B2) becomes 



(B3) 


For the logarithmic decrement, equation (B3) at t = 0 becomes 


a |o = | a 


and, after n cycles, at t = 2jm 


i/o>i. 


a| n = a|e 


-TtAn 


or 


l a |o 


a 


= e 


JtAn 


n 


(b4) 


By taking the natural logarithm of both sides of equation (B4), 




" Bt 


which is measured with the wing subjected to air flow. Returning to 
the definition of A, 


k s g s - | pv 2 c 2 sm 2 
St = 

k g - ^ pv 2 c 2 sm 1 



(B5) 
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Since 


k s ~ | pv 2 c 2 sm 1 = Igo^ 2 


and 


k s = I s«Vac 2 


then 


-m2 = 


I S 0 > I 2 


2 


p-^c^s 




(b6) 


Equation (b 6) is the form used in the reduction of the damping -moment 
data in this paper. 
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APPENDIX C 

SOME EXPERIMENTAL EVIDENCE OF WIND -TUNNEL-WALL EFFECTS 
ON AN OSCILLATING WING 


This appendix is presented to document data relating to some wind- 
tunnel-wall effects on oscillating air forces and moments. The problem 
of the effect of the presence of wind-tunnel walls on measurements per- 
taining to airfoils and wings in the steady-state case has_ been resolved. 
The problem for the condition of unsteady flow has been treated theo- 
retically by Reissner in reference 11, Jones in reference 12, and Timman 
in reference 13 for the incompressible case, and by Runyan and Watkins 
in reference 9 for the two-dimensional subsonic compressible case. For 
the incompressible -flow condition, the influence of the tunnel walls has 
been found to be comparatively small for most cases , although indications 
are given that, for some ranges of 2v/ca>, the effect may be very large. 
For the compressible-flow case, reference 14- indicates the possibility 
of obtaining a resonant condition which might result in a misinterpre- 
tation of the measured quantities, the critical condition for the rec- 
tangular tunnel being given by oil /a = (2m - l.)jt \Jl - where m = 1, 

2, and 3* 

Reference 9 indicates that this resonant condition corresponds to 
the establishment of transverse velocities having a maximum amplitude 
at the airfoil. These transverse velocity components alter the effec- 
tive angle of attack and thus affect the air forces. Such a condition 
might be expected to be obtained for the finite wing, although the 
behavior of the measured forces and moments in approaching the vicinity 
of the "critical" condition may not necessarily follow the pattern of 
the two-dimensional wing. 

An indication that a distortion and a resonance is experienced by 
oscillating wings has been obtained experimentally. These data are pre- 
sented in table VIII and figures 10 and 11. It might be remarked that 
the phenomena have also been observed for the two-dimensional wing in 
connection with the work of reference 5. The peculiar behavior of the 
data in the neighborhood of the Mach number critical for the resonance 
condition is strikingly demonstrated by the damping-moment coefficients 
(see fig. 10). For the finite wing, these data tend to yield peak values 
of the damping-moment coefficients in the vicinity of the "critical Mach 
number." It is of interest to note that this behavior is the inverse 
to that observed for the data of the infinite wing of reference 13 . 

This inversion may be attributed to the introduction of three-dimensional 
effects into the problem. 
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Figure 11 indicates effects of wall interference on the oscillating 
lift coefficients. As the critical value of l/k is reached, the lift 
coefficient drops off and there is a definite dip in the curve. 

These data have been presented to furnish some evidence of the 
nature of the effects of tunnel walls on the aerodynamic forces and 
moments acting on an oscillating wing. It is evident that some caution 
is required to avoid conditions leading to these effects. 
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TABLE II.- EXPERIMENTAL DATA FOR TORSIONAL SPRING B 
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TABLE V.- THEORETICAL INC(MPRESSIBLE COEFFICIENTS OF 


REISSNER APPLIED AT ASPECT RATIO 2 



|L| 


jtpv^sc a 



6.67 8.08 

5.00 11.01 

4.00 14.50 
2.50 24.20 

1.67 33.90 
1.25 44.70 

1.00 52.50 



0.9320 O.1315 - 10.25 0.4684 

.9344 .1819 - 9.15 .4704 

.9244 .2396 -15.28 .4736 

.9288 .3984 -21.80 .5044 

.9530 .6400 -26.00 .5644 

.9272 .9184 -28.70 .6232 

.9564 1.2456 -27.30 .7292 


TABLE VI.- THEORETICAL INCOMPRESSIBLE COEFFICIENTS OF 
LAWRENCE AND GERBER FOR ASPECT RATIO 2 
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TABLE VII.- TWO-DIMENSIONAL INCOMPRESSIBLE COEFFICIENTS OF THEODORSEN 
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TABLE VTII. - EXPERIMENTAL DATA SHOWING EFFECT OF WIND-TUNNEL-WALL 
INTERFERENCE FOR TORSIONAL SPRING D' IN FREON 
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Figure 1.- Diagram of oscillating mechanism and wing mounted in the 
Langley 2- by 4-foot flutter research tunnel. 
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Figure 2.- Experimental aerodynamic coefficients against k, showing 
comparison of data of this report with data of reference 6. For 

reference 6 data, M = 0.13 and 0.7 X 10^ < R < 0.88 x 10^. 
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Figure 4.- Theoretical and experimental lift coefficients against l/k. 
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Figure 6.- Theoretical and experimental inphase-moment-component 
coefficients against l/k for torsional springs A and B. 
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Figure 9 • - Reduced frequency against Mach number showing range of 

experimental studies. 
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Figure 11.- Theoretical and experimental lift coefficients indicating 
tunnel-wall effects plotted against l/k for torsion spring D' . 
M«0.35(l/k). 





